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CONSTRUCTION OF THREE DIMENSIONAL SOLUTIONS FOR THE MAXWELL 

EQUATIONS* 

A. YEFET* AND E. TURKEL* 

Abstract. Wc consider numerical solutions for the three dimensional time dependent Maxwell equations. 
We construct a fourth order accurate compact implicit scheme and compare it to the Ycc scheme for free 
space in a box. 

Subject classification. Applied and Numerical Mathematics 
Key words. Maxwell equations, the Yce scheme, the Ty(2,4) scheme 

1. Maxwell Equations in a Box. Let r = ci = tj and Z — For the rest of this paper wc 

replace r by t. The three dimensional time dependent Maxwell equations then arc: 

dE x =z (dH ± _ dH y \ 
dt \ dy dz ) 

dE y =z (dH x dH z \ 
dt \ dz dx ) 

dE z =z (dHy dH A 

dt \ dx dy ) 

dH x 1 (dE y dE z \ 

~df ~ Z ( dz ~ dy ) 
dHy _ 1 ( dE z _ dE x \ 
dt Z \ dx dz ) 
dHz^l/dE* dE y \ 

[ } dt Z\dy dx )' 

We set Z = 1 in this paper. 

A plane wave solution is given by 

H x = sin(wt) sin( Ax + Dy + Cz) 

H y = Hy sin (ujt) sin(Ax 4- By + Cz) 

H z = H® sin(o;t) sin(Ar -f By -h Cz) 

E x = E x cos(u ft) cos(Ar -j- By + Cz) 

E y = Ey cos(u jt) cos (Ax + By -j- Cz) 

E z = E® cos (ut) cos( Ax -f By -F Cz) 

Substituting into the Maxwell equations this is a solution if 

(1.2) uj 2 = A 2 + B 2 + C 2 

(1.3) 0 = AH* + BH° + CH° Z 
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We also demand that 


loE° x = HyC - H° Z B 
u>E° = H° Z A - H° X C 
ujE° z = H°B - H°A 

2. Numerical Tests. We consider a case where H ? = H?. = = 1 and 

x y * 


A = 7T 
B = —2n 
C — 7T 

U) = V67T 


We use this exact solution as a basis for comparison in the box [0, 1/2] x [0, 1/4] x [0, 1/2] . We shall compare 
two numerical methods: the Yec scheme [1] which is second order accurate in space and time and the Ty(2,4) 
scheme [2, 3] which is second order accurate in time but fourth order accurate in space. In order for the 
total error to be fourth order we must choose the time step small enough so that the temporal error does 
not swamp the spatial error. This requires At ~ (Ax) 2 . If the error requirements are too severe then this 
is inefficient and the leapfrog in time should be replaced by a fourth order Runge-Kutta method. However, 
for the experiments in this paper we shall use the same leapfrog method for both schemes. Hence, both the 
Yec scheme and the Ty(2,4) have the electric and magnetic variables at the same staggered locations both 
in space and in time. The Yec scheme approximates the derivatives via the following approximation. 
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A similar formula holds for the other variables shifted to other locations in each direction. The Ty(2,4) 
scheme is an implicit compact scheme given by 
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where A is defined the following way: 
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For the Yee scheme we choose At = ^ while for thc T y( 2 > 4 ) scheme we choose At ~ h 2 where h = 
Ax ~ Ay. 

We measure the error in thc L 2 norm between the approximate and exact electric field in thc z-direction. 
Thc Ty(2,4) behaves better than expected and gives almost fifth order accuracy. The Yee scheme gives a 
second order accuracy as expected. 


REFERENCES 

[1] K.S. YEE, Numerical Solution of Initial Boundary Value Problems Involving Maxwell's Equation in 

Isotropic Media , IEEE Transactions on Antennas and Propagation AP-14 (May 1996), pp. 302- 
307. 

[2] E. Turkel and A. Yefet, Fourth Order Accurate Compact Implicit Method for the Maxwell Equa- 

tions , submitted to IEEE Transactions on Antennas and Propagation. 

[3] A. Yefet and E. Turkel, Fourth Order Compact Implicit Method for the Maxwell Equations with 

Discontinuous Coefficients, submitted to Applied Numerical Mathematics. 


3 



scheme 

h 

At 

t= 10 

reduction 
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Ty{ 2,4) 
Ty( 2,4) 
Ty( 2,4) 
Fee 
Yee 
Yee 

1 

20 

1 

40 

1 

80 

1 

20 

1 

40 

1 

80 

1 
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1 

1600 

1 
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1 

35 

1 

70 

1 

140 

3.62 x 1(T 4 
1.1443 x 10“ 5 
3.5621 x 1CT 7 
0.027 
7.3 x 10“ 4 
1.82 x 10^ 4 

31.6423 

32.1255 

3.694 

4.0042 

4.98 

5.0056 

1.9028 
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Table 2.1 

Comparison of the maximum errors in L2 norm 



Fig. 2.1. log! 0 ( error) For the Yee scheme. 
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